3.1.57 \(\int \frac {x^5}{(a x^2+b x^3+c x^4)^{3/2}} \, dx\) [57]

3.1.57.1 Optimal result
3.1.57.2 Mathematica [A] (verified)
3.1.57.3 Rubi [A] (verified)
3.1.57.4 Maple [A] (verified)
3.1.57.5 Fricas [A] (verification not implemented)
3.1.57.6 Sympy [F]
3.1.57.7 Maxima [F]
3.1.57.8 Giac [A] (verification not implemented)
3.1.57.9 Mupad [F(-1)]

3.1.57.1 Optimal result

Integrand size = 24, antiderivative size = 153 \[ \int \frac {x^5}{\left (a x^2+b x^3+c x^4\right )^{3/2}} \, dx=\frac {2 x^2 (2 a+b x)}{\left (b^2-4 a c\right ) \sqrt {a x^2+b x^3+c x^4}}-\frac {2 b \sqrt {a x^2+b x^3+c x^4}}{c \left (b^2-4 a c\right ) x}+\frac {x \sqrt {a+b x+c x^2} \text {arctanh}\left (\frac {b+2 c x}{2 \sqrt {c} \sqrt {a+b x+c x^2}}\right )}{c^{3/2} \sqrt {a x^2+b x^3+c x^4}} \]

output
2*x^2*(b*x+2*a)/(-4*a*c+b^2)/(c*x^4+b*x^3+a*x^2)^(1/2)+x*arctanh(1/2*(2*c* 
x+b)/c^(1/2)/(c*x^2+b*x+a)^(1/2))*(c*x^2+b*x+a)^(1/2)/c^(3/2)/(c*x^4+b*x^3 
+a*x^2)^(1/2)-2*b*(c*x^4+b*x^3+a*x^2)^(1/2)/c/(-4*a*c+b^2)/x
 
3.1.57.2 Mathematica [A] (verified)

Time = 0.35 (sec) , antiderivative size = 112, normalized size of antiderivative = 0.73 \[ \int \frac {x^5}{\left (a x^2+b x^3+c x^4\right )^{3/2}} \, dx=-\frac {x \left (2 \sqrt {c} \left (-a b-b^2 x+2 a c x\right )+\left (b^2-4 a c\right ) \sqrt {a+x (b+c x)} \text {arctanh}\left (\frac {b+2 c x}{2 \sqrt {c} \sqrt {a+x (b+c x)}}\right )\right )}{c^{3/2} \left (-b^2+4 a c\right ) \sqrt {x^2 (a+x (b+c x))}} \]

input
Integrate[x^5/(a*x^2 + b*x^3 + c*x^4)^(3/2),x]
 
output
-((x*(2*Sqrt[c]*(-(a*b) - b^2*x + 2*a*c*x) + (b^2 - 4*a*c)*Sqrt[a + x*(b + 
 c*x)]*ArcTanh[(b + 2*c*x)/(2*Sqrt[c]*Sqrt[a + x*(b + c*x)])]))/(c^(3/2)*( 
-b^2 + 4*a*c)*Sqrt[x^2*(a + x*(b + c*x))]))
 
3.1.57.3 Rubi [A] (verified)

Time = 0.34 (sec) , antiderivative size = 166, normalized size of antiderivative = 1.08, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.250, Rules used = {1970, 1996, 27, 1961, 1092, 219}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {x^5}{\left (a x^2+b x^3+c x^4\right )^{3/2}} \, dx\)

\(\Big \downarrow \) 1970

\(\displaystyle \frac {2 x^2 (2 a+b x)}{\left (b^2-4 a c\right ) \sqrt {a x^2+b x^3+c x^4}}-\frac {2 \int \frac {x (2 a+b x)}{\sqrt {c x^4+b x^3+a x^2}}dx}{b^2-4 a c}\)

\(\Big \downarrow \) 1996

\(\displaystyle \frac {2 x^2 (2 a+b x)}{\left (b^2-4 a c\right ) \sqrt {a x^2+b x^3+c x^4}}-\frac {2 \left (\frac {b \sqrt {a x^2+b x^3+c x^4}}{c x}-\frac {\int \frac {\left (b^2-4 a c\right ) x}{2 \sqrt {c x^4+b x^3+a x^2}}dx}{c}\right )}{b^2-4 a c}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {2 x^2 (2 a+b x)}{\left (b^2-4 a c\right ) \sqrt {a x^2+b x^3+c x^4}}-\frac {2 \left (\frac {b \sqrt {a x^2+b x^3+c x^4}}{c x}-\frac {\left (b^2-4 a c\right ) \int \frac {x}{\sqrt {c x^4+b x^3+a x^2}}dx}{2 c}\right )}{b^2-4 a c}\)

\(\Big \downarrow \) 1961

\(\displaystyle \frac {2 x^2 (2 a+b x)}{\left (b^2-4 a c\right ) \sqrt {a x^2+b x^3+c x^4}}-\frac {2 \left (\frac {b \sqrt {a x^2+b x^3+c x^4}}{c x}-\frac {x \left (b^2-4 a c\right ) \sqrt {a+b x+c x^2} \int \frac {1}{\sqrt {c x^2+b x+a}}dx}{2 c \sqrt {a x^2+b x^3+c x^4}}\right )}{b^2-4 a c}\)

\(\Big \downarrow \) 1092

\(\displaystyle \frac {2 x^2 (2 a+b x)}{\left (b^2-4 a c\right ) \sqrt {a x^2+b x^3+c x^4}}-\frac {2 \left (\frac {b \sqrt {a x^2+b x^3+c x^4}}{c x}-\frac {x \left (b^2-4 a c\right ) \sqrt {a+b x+c x^2} \int \frac {1}{4 c-\frac {(b+2 c x)^2}{c x^2+b x+a}}d\frac {b+2 c x}{\sqrt {c x^2+b x+a}}}{c \sqrt {a x^2+b x^3+c x^4}}\right )}{b^2-4 a c}\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {2 x^2 (2 a+b x)}{\left (b^2-4 a c\right ) \sqrt {a x^2+b x^3+c x^4}}-\frac {2 \left (\frac {b \sqrt {a x^2+b x^3+c x^4}}{c x}-\frac {x \left (b^2-4 a c\right ) \sqrt {a+b x+c x^2} \text {arctanh}\left (\frac {b+2 c x}{2 \sqrt {c} \sqrt {a+b x+c x^2}}\right )}{2 c^{3/2} \sqrt {a x^2+b x^3+c x^4}}\right )}{b^2-4 a c}\)

input
Int[x^5/(a*x^2 + b*x^3 + c*x^4)^(3/2),x]
 
output
(2*x^2*(2*a + b*x))/((b^2 - 4*a*c)*Sqrt[a*x^2 + b*x^3 + c*x^4]) - (2*((b*S 
qrt[a*x^2 + b*x^3 + c*x^4])/(c*x) - ((b^2 - 4*a*c)*x*Sqrt[a + b*x + c*x^2] 
*ArcTanh[(b + 2*c*x)/(2*Sqrt[c]*Sqrt[a + b*x + c*x^2])])/(2*c^(3/2)*Sqrt[a 
*x^2 + b*x^3 + c*x^4])))/(b^2 - 4*a*c)
 

3.1.57.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 219
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))* 
ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && NegQ[a/b] && (Gt 
Q[a, 0] || LtQ[b, 0])
 

rule 1092
Int[1/Sqrt[(a_) + (b_.)*(x_) + (c_.)*(x_)^2], x_Symbol] :> Simp[2   Subst[I 
nt[1/(4*c - x^2), x], x, (b + 2*c*x)/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a 
, b, c}, x]
 

rule 1961
Int[(x_)^(m_.)/Sqrt[(b_.)*(x_)^(n_.) + (a_.)*(x_)^(q_.) + (c_.)*(x_)^(r_.)] 
, x_Symbol] :> Simp[x^(q/2)*(Sqrt[a + b*x^(n - q) + c*x^(2*(n - q))]/Sqrt[a 
*x^q + b*x^n + c*x^(2*n - q)])   Int[x^(m - q/2)/Sqrt[a + b*x^(n - q) + c*x 
^(2*(n - q))], x], x] /; FreeQ[{a, b, c, m, n, q}, x] && EqQ[r, 2*n - q] && 
 PosQ[n - q] && ((EqQ[m, 1] && EqQ[n, 3] && EqQ[q, 2]) || ((EqQ[m + 1/2] || 
 EqQ[m, 3/2] || EqQ[m, 1/2] || EqQ[m, 5/2]) && EqQ[n, 3] && EqQ[q, 1]))
 

rule 1970
Int[(x_)^(m_.)*((b_.)*(x_)^(n_.) + (a_.)*(x_)^(q_.) + (c_.)*(x_)^(r_.))^(p_ 
), x_Symbol] :> Simp[(-x^(m - 2*n + q + 1))*(2*a + b*x^(n - q))*((a*x^q + b 
*x^n + c*x^(2*n - q))^(p + 1)/((n - q)*(p + 1)*(b^2 - 4*a*c))), x] + Simp[1 
/((n - q)*(p + 1)*(b^2 - 4*a*c))   Int[x^(m - 2*n + q)*(2*a*(m + p*q - 2*(n 
 - q) + 1) + b*(m + p*q + (n - q)*(2*p + 1) + 1)*x^(n - q))*(a*x^q + b*x^n 
+ c*x^(2*n - q))^(p + 1), x], x] /; FreeQ[{a, b, c}, x] && EqQ[r, 2*n - q] 
&& PosQ[n - q] &&  !IntegerQ[p] && NeQ[b^2 - 4*a*c, 0] && IGtQ[n, 0] && LtQ 
[p, -1] && RationalQ[m, q] && GtQ[m + p*q + 1, 2*(n - q)]
 

rule 1996
Int[(x_)^(m_.)*((c_.)*(x_)^(j_.) + (b_.)*(x_)^(n_.) + (a_.)*(x_)^(q_.))^(p_ 
.)*((A_) + (B_.)*(x_)^(r_.)), x_Symbol] :> Simp[B*x^(m - n + 1)*((a*x^q + b 
*x^n + c*x^(2*n - q))^(p + 1)/(c*(m + p*q + (n - q)*(2*p + 1) + 1))), x] - 
Simp[1/(c*(m + p*q + (n - q)*(2*p + 1) + 1))   Int[x^(m - n + q)*Simp[a*B*( 
m + p*q - n + q + 1) + (b*B*(m + p*q + (n - q)*p + 1) - A*c*(m + p*q + (n - 
 q)*(2*p + 1) + 1))*x^(n - q), x]*(a*x^q + b*x^n + c*x^(2*n - q))^p, x], x] 
 /; FreeQ[{a, b, c, A, B}, x] && EqQ[r, n - q] && EqQ[j, 2*n - q] &&  !Inte 
gerQ[p] && NeQ[b^2 - 4*a*c, 0] && IGtQ[n, 0] && GeQ[p, -1] && LtQ[p, 0] && 
RationalQ[m, q] && GeQ[m + p*q, n - q - 1] && NeQ[m + p*q + (n - q)*(2*p + 
1) + 1, 0]
 
3.1.57.4 Maple [A] (verified)

Time = 0.15 (sec) , antiderivative size = 84, normalized size of antiderivative = 0.55

method result size
pseudoelliptic \(-\frac {x}{c \sqrt {c \,x^{2}+b x +a}}+\frac {b \left (b x +2 a \right )}{\sqrt {c \,x^{2}+b x +a}\, c \left (4 a c -b^{2}\right )}+\frac {\ln \left (2 \sqrt {c \,x^{2}+b x +a}\, \sqrt {c}+2 c x +b \right )}{c^{\frac {3}{2}}}\) \(84\)
default \(-\frac {x^{3} \left (c \,x^{2}+b x +a \right ) \left (4 c^{\frac {5}{2}} a x -2 c^{\frac {3}{2}} b^{2} x -2 c^{\frac {3}{2}} a b -4 \ln \left (\frac {2 \sqrt {c \,x^{2}+b x +a}\, \sqrt {c}+2 c x +b}{2 \sqrt {c}}\right ) \sqrt {c \,x^{2}+b x +a}\, a \,c^{2}+\ln \left (\frac {2 \sqrt {c \,x^{2}+b x +a}\, \sqrt {c}+2 c x +b}{2 \sqrt {c}}\right ) \sqrt {c \,x^{2}+b x +a}\, b^{2} c \right )}{c^{\frac {5}{2}} \left (c \,x^{4}+b \,x^{3}+a \,x^{2}\right )^{\frac {3}{2}} \left (4 a c -b^{2}\right )}\) \(166\)

input
int(x^5/(c*x^4+b*x^3+a*x^2)^(3/2),x,method=_RETURNVERBOSE)
 
output
-x/c/(c*x^2+b*x+a)^(1/2)+b*(b*x+2*a)/(c*x^2+b*x+a)^(1/2)/c/(4*a*c-b^2)+1/c 
^(3/2)*ln(2*(c*x^2+b*x+a)^(1/2)*c^(1/2)+2*c*x+b)
 
3.1.57.5 Fricas [A] (verification not implemented)

Time = 0.30 (sec) , antiderivative size = 414, normalized size of antiderivative = 2.71 \[ \int \frac {x^5}{\left (a x^2+b x^3+c x^4\right )^{3/2}} \, dx=\left [\frac {{\left ({\left (b^{2} c - 4 \, a c^{2}\right )} x^{3} + {\left (b^{3} - 4 \, a b c\right )} x^{2} + {\left (a b^{2} - 4 \, a^{2} c\right )} x\right )} \sqrt {c} \log \left (-\frac {8 \, c^{2} x^{3} + 8 \, b c x^{2} + 4 \, \sqrt {c x^{4} + b x^{3} + a x^{2}} {\left (2 \, c x + b\right )} \sqrt {c} + {\left (b^{2} + 4 \, a c\right )} x}{x}\right ) - 4 \, \sqrt {c x^{4} + b x^{3} + a x^{2}} {\left (a b c + {\left (b^{2} c - 2 \, a c^{2}\right )} x\right )}}{2 \, {\left ({\left (b^{2} c^{3} - 4 \, a c^{4}\right )} x^{3} + {\left (b^{3} c^{2} - 4 \, a b c^{3}\right )} x^{2} + {\left (a b^{2} c^{2} - 4 \, a^{2} c^{3}\right )} x\right )}}, -\frac {{\left ({\left (b^{2} c - 4 \, a c^{2}\right )} x^{3} + {\left (b^{3} - 4 \, a b c\right )} x^{2} + {\left (a b^{2} - 4 \, a^{2} c\right )} x\right )} \sqrt {-c} \arctan \left (\frac {\sqrt {c x^{4} + b x^{3} + a x^{2}} {\left (2 \, c x + b\right )} \sqrt {-c}}{2 \, {\left (c^{2} x^{3} + b c x^{2} + a c x\right )}}\right ) + 2 \, \sqrt {c x^{4} + b x^{3} + a x^{2}} {\left (a b c + {\left (b^{2} c - 2 \, a c^{2}\right )} x\right )}}{{\left (b^{2} c^{3} - 4 \, a c^{4}\right )} x^{3} + {\left (b^{3} c^{2} - 4 \, a b c^{3}\right )} x^{2} + {\left (a b^{2} c^{2} - 4 \, a^{2} c^{3}\right )} x}\right ] \]

input
integrate(x^5/(c*x^4+b*x^3+a*x^2)^(3/2),x, algorithm="fricas")
 
output
[1/2*(((b^2*c - 4*a*c^2)*x^3 + (b^3 - 4*a*b*c)*x^2 + (a*b^2 - 4*a^2*c)*x)* 
sqrt(c)*log(-(8*c^2*x^3 + 8*b*c*x^2 + 4*sqrt(c*x^4 + b*x^3 + a*x^2)*(2*c*x 
 + b)*sqrt(c) + (b^2 + 4*a*c)*x)/x) - 4*sqrt(c*x^4 + b*x^3 + a*x^2)*(a*b*c 
 + (b^2*c - 2*a*c^2)*x))/((b^2*c^3 - 4*a*c^4)*x^3 + (b^3*c^2 - 4*a*b*c^3)* 
x^2 + (a*b^2*c^2 - 4*a^2*c^3)*x), -(((b^2*c - 4*a*c^2)*x^3 + (b^3 - 4*a*b* 
c)*x^2 + (a*b^2 - 4*a^2*c)*x)*sqrt(-c)*arctan(1/2*sqrt(c*x^4 + b*x^3 + a*x 
^2)*(2*c*x + b)*sqrt(-c)/(c^2*x^3 + b*c*x^2 + a*c*x)) + 2*sqrt(c*x^4 + b*x 
^3 + a*x^2)*(a*b*c + (b^2*c - 2*a*c^2)*x))/((b^2*c^3 - 4*a*c^4)*x^3 + (b^3 
*c^2 - 4*a*b*c^3)*x^2 + (a*b^2*c^2 - 4*a^2*c^3)*x)]
 
3.1.57.6 Sympy [F]

\[ \int \frac {x^5}{\left (a x^2+b x^3+c x^4\right )^{3/2}} \, dx=\int \frac {x^{5}}{\left (x^{2} \left (a + b x + c x^{2}\right )\right )^{\frac {3}{2}}}\, dx \]

input
integrate(x**5/(c*x**4+b*x**3+a*x**2)**(3/2),x)
 
output
Integral(x**5/(x**2*(a + b*x + c*x**2))**(3/2), x)
 
3.1.57.7 Maxima [F]

\[ \int \frac {x^5}{\left (a x^2+b x^3+c x^4\right )^{3/2}} \, dx=\int { \frac {x^{5}}{{\left (c x^{4} + b x^{3} + a x^{2}\right )}^{\frac {3}{2}}} \,d x } \]

input
integrate(x^5/(c*x^4+b*x^3+a*x^2)^(3/2),x, algorithm="maxima")
 
output
integrate(x^5/(c*x^4 + b*x^3 + a*x^2)^(3/2), x)
 
3.1.57.8 Giac [A] (verification not implemented)

Time = 0.35 (sec) , antiderivative size = 164, normalized size of antiderivative = 1.07 \[ \int \frac {x^5}{\left (a x^2+b x^3+c x^4\right )^{3/2}} \, dx=\frac {{\left (b^{2} \log \left ({\left | b - 2 \, \sqrt {a} \sqrt {c} \right |}\right ) - 4 \, a c \log \left ({\left | b - 2 \, \sqrt {a} \sqrt {c} \right |}\right ) + 2 \, \sqrt {a} b \sqrt {c}\right )} \mathrm {sgn}\left (x\right )}{b^{2} c^{\frac {3}{2}} - 4 \, a c^{\frac {5}{2}}} - \frac {2 \, {\left (\frac {a b}{b^{2} c \mathrm {sgn}\left (x\right ) - 4 \, a c^{2} \mathrm {sgn}\left (x\right )} + \frac {{\left (b^{2} - 2 \, a c\right )} x}{b^{2} c \mathrm {sgn}\left (x\right ) - 4 \, a c^{2} \mathrm {sgn}\left (x\right )}\right )}}{\sqrt {c x^{2} + b x + a}} - \frac {\log \left ({\left | 2 \, {\left (\sqrt {c} x - \sqrt {c x^{2} + b x + a}\right )} \sqrt {c} + b \right |}\right )}{c^{\frac {3}{2}} \mathrm {sgn}\left (x\right )} \]

input
integrate(x^5/(c*x^4+b*x^3+a*x^2)^(3/2),x, algorithm="giac")
 
output
(b^2*log(abs(b - 2*sqrt(a)*sqrt(c))) - 4*a*c*log(abs(b - 2*sqrt(a)*sqrt(c) 
)) + 2*sqrt(a)*b*sqrt(c))*sgn(x)/(b^2*c^(3/2) - 4*a*c^(5/2)) - 2*(a*b/(b^2 
*c*sgn(x) - 4*a*c^2*sgn(x)) + (b^2 - 2*a*c)*x/(b^2*c*sgn(x) - 4*a*c^2*sgn( 
x)))/sqrt(c*x^2 + b*x + a) - log(abs(2*(sqrt(c)*x - sqrt(c*x^2 + b*x + a)) 
*sqrt(c) + b))/(c^(3/2)*sgn(x))
 
3.1.57.9 Mupad [F(-1)]

Timed out. \[ \int \frac {x^5}{\left (a x^2+b x^3+c x^4\right )^{3/2}} \, dx=\int \frac {x^5}{{\left (c\,x^4+b\,x^3+a\,x^2\right )}^{3/2}} \,d x \]

input
int(x^5/(a*x^2 + b*x^3 + c*x^4)^(3/2),x)
 
output
int(x^5/(a*x^2 + b*x^3 + c*x^4)^(3/2), x)